Harmonic oscillator in a rotating trap: Complete solution in 3D 
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Complete description of the classical and quantum dynamics of a particle in an anisotropic, 
rotating, harmonic trap is given. The problem is studied in three dimensions and no restrictions on 
the geometry are imposed. In the generic case, for an arbitrary orientation of the rotation axis, there 
are two regions of instability with different characteristics. The analysis of the quantum-mechanical 
problem is made simple due to a direct connection between the classical mode vectors and the 
quantum-mechanical wave functions. This connection is obtained via the matrix Riccati equation 
that governs the time evolution of squeezed states of the harmonic oscillator. It is also shown that 
the inclusion of gravity leads to a resonant behavior — the particles are expelled from the trap. 
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I. INTRODUCTION 

Harmonic traps are often used in optics and atomic 
physics (especially in the study of Bose-Einstein conden- 
sates) and yet a complete theory of these devices has not 
been developed. In this work we present a complete so- 
lution to the problem of the motion of a single particle 
classical or quantum-mechanical moving in three dimen- 
sions for the most general anisotropic rotating harmonic 
trap. This is an exactly soluble problem but technical 
difficulties apparently served so far as a deterrent in de- 
veloping a full description. All our results are valid also 
for the center of mass motion in many-body theory since 
for the quadratic potentials the center of mass motion 
completely separates for the internal motion 0, 0, . 

The solution to the problem of a two-dimensional ro- 
tating trap has been known for at least hundred years. 
In the classic textbook on analytical dynamics by Whit- 
taker 01 one finds a solution of a mathematically equiv- 
alent problem of small oscillations of "a heavy particle 
about its position of equilibrium at the lowest point of a 
surface which is rotating with constant angular velocity 
about a vertical axis through the point" . Recently, this 
problem acquired a new significance in connection with 
the experimental and theoretical studies of Bose-Einstein 
condensates in rotating traps 0, IE IE S ■ Still, to our 
knowledge, a full description of the particle dynamics in 
a rotating anisotropic trap in three dimensions was not 
given. Explicit formulas describing the complete mode 
structure in the three-dimensional case are indeed quite 
cumbersome |Tol | because we deal here with third-order 
polynomials and on top of that they have rather com- 
plicated coefficients. However, many important features 
may be exhibited without straining the reader's patience. 
In particular, we determine the characteristic frequencies 
for an arbitrary orientation of the angular velocity, and 



we identify various stability regions. Following a recent 
theoretical observation that the gravity may cause reso- 
nances |llj , we give a full analysis of these phenomenon. 
Finally, owing to the properties of the matrix Riccati 
equations we unravel a direct connection between the so- 
lutions of the classical problem and the shape of the quan- 
tum wave functions. Also, we find the three constants of 
motion quadratic in the positions and momenta. 



II. EQUATIONS OF MOTION 

The solution of our problem is the simplest in the co- 
rotating coordinate system, where the potential does not 
depend on time and the Hamiltonian has the following 
form: 

- m - 
li = +rVlp+ —rVr. 1) 
2m 2 

The antisymmetric matrix Cl is built from the compo- 
nents of the angular velocity vector 17^^ = Cijk^j and the 
potential matrix V is symmetric and positive definite. 
The term involving is responsible for both the Coriolis 
force and the centrifugal force. In quantum mechanics 
the Hamiltonian ^ may be obtained by transforming 
the Schrodinger equation to the rotating frame by a sub- 
stitution ij} — > exp(in-Lt/h)'ip. 

The equations of motion determined by look the 
same in classical and in quantum theory and they can be 
written in the following matrix form 

^ = MX(i), (2) 
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where the six-dimensional vector X(t) comprises coordi- 
nates and momenta X = (cc, y, z,px,Py,Pz) and the 6x6 
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matrix M is 
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The set of six linear differential equations Q has in gen- 
eral six independent solutions Xi{t) (normal modes) that 
in the stability region change harmonically with time. 
The frequencies of the normal modes are determined by 
the characteristic polynomial P{uj) associated with the 
matrix M. Owing to the time-reversal symmetry of our 
problem, the characteristic polynomial is triquadratic 
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FIG. 1: The contour lines in this figure represent the zeroes of 
the characteristic polynomial JSJ for Vx — 1, Vy — 2, Vz — 3 
and Tlx = l/\/3, riy = l/VS, = l/VS plotted as functions 



The coefficients A, B, C have been obtained in our previ- 
ous publication jl^] and have the following rotationally 
invariant forms 



the magnitude of angular velocity Q and the square x ■ 



of 



A = -2fl^-Ti{V}, 

B ^ fl'^ + fl^{3n-V-n-TT{V})^ 



Tr{y}2-Tr{l>2} 



(5a) 
,(5b) 



the characteristic frequency. In addition to the lower region 
of instability (between dashed lines) where one of the roots of 
@ is negative, there is also the upper region (between dotted 
lines) where there exists only one real root (described by the 
continuation of the line that begins at x ~ 3, not seen in this 
figure because it corresponds to a very large value of x)- 



C = fl''{Tv{V}-fl^)n-V-n-fl^n-V^-n-Dct{V},{5c) expressions in 



where fl is the length and n is the direction of the an- 
gular velocity vector, fl = fin. Since the matrix V 
is by definition symmetric, one may always choose a 
base in the rotating frame such that V will be diagonal 
{V = diag(14, V^, Vz))- As was to be expected, the mass 
has dropped out completely from the characteristic equa- 
tion. Therefore, all the results will scale uniformly and 
one may express them in terms of an arbitrarily chosen 
frequency unit t^o- In what follows, all the dimensionless 
frequencies are to be the multiplied by ujq and the trap 
parameters Vi by Wq. 



III. STABILITY OF MOTION 

The stability of motion for a harmonic oscillator is de- 
termined by the values of its characteristic frequencies 
w — the roots of the characteristic polynomial. Stable 
oscillations take place when all w's are real. This means 
that all three roots of the polynomial Q{x) 



fl± = 



'b±Vb'^ - 4ac 



2a 



(7) 



Qix)^X^ + Ax^ + Bx + C, x = c^^ 



(6) 



are real and positive. It has been argued in Ref. ^Td\ that 
there is always a region of instability; this occurs when 
one of the roots of Q{x) is negative. The corresponding 
modes grow exponentially with time. As seen in Fig. ^ 
this region of instability is bounded by the two values fl± 
at which the curve crosses the vertical axis. These values 
are given by the zeroes of C, treated as a biquadratic 



where a = nV n, b = Tr{V}n -V -n — nV^ -n, and 
c = Det{F}. Since a, b, and c are positive, fl± are real. 

In this paper we show that a more detailed investiga- 
tion reveals that in addition to the region of (purely ex- 
ponential) instability when one root of Q{x) is negative 
there is also, in general, an additional region of (oscilla- 
tory) instability where two roots are complex. Without 
rotation, when = 0, the three roots of Q{x) are equal 
to the eigenvalues of the potential matrix V and one has 
a simple system of three harmonic oscillators vibrating 
independently along the principal directions of the trap. 
As fl increases, our system has, in general, two regions of 
instability: the lower one when one of the roots of Q{x) 
is negative and the upper one when two roots are com- 
plex. We shall exhibit this behavior by plotting the zero 
contour lines of lO in the flx-p^a^ne. We assume that the 
trap potential and the direction of rotation are fixed and 
we treat the characteristic polynomial Q{x) a-s a function 
of fl and X only. Contour lines representing the zeroes 
of JBJ in the generic case are shown in Fig. ^ There is 
a region of fl, where only one real solution of Q exists. 
However, this region is bounded, so for sufficiently large 
fl the system is always stable. 

Graphical representation of the solutions for Ux — 1, 
Uy — = is shown in Fig. |21 In this plot only one 
instability region is seen, where one of the roots of © is 
negative. Owing to the stabilizing effect of the Coriolis 
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FIG. 2: This plot is for the same trap as in Fig. but for 
the rotation around the trap axis Ux = 0, Uy = 0, Uz — 1. In 
this degenerate case there exists only one region of instability, 
when one of the roots of ijSJ is negative. 
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FIG. 3: This plot shows how the degenerate case (Fig. I^J 
merges with the general case (Fig. [Hi. Here we have chosen 
the following values of the parameters: Vx = 1, Vy = 2, Vz ~ 3 
and rix = sin(l/10), Uy = 0, Uz = cos(l/10). 



force, for fast rotations the system becomes again stable. 
In this plot one can also see, that one of the character- 
istic frequencies remains constant. It is so, because the 
rotation does not influence the motion in the direction 
of the rotation axis. When fl increases, the solutions of 
© change their character. In the previous case there 
was only one instability region, where the square of the 
frequency was negative. When the direction of angular 
velocity is not parallel to one of the axes of the trap, a 
second kind of instability appears — the square of the 
frequency becomes complex. Since the coefficients of the 
characteristic polynomial are real, the square of the sec- 
ond frequency is complex conjugate to the first one. In 
this case the instability has the form of expanding oscilla- 
tions. In Fig.|21we show what happens when the rotation 
axis is tilted away a little bit from the z direction and the 
second kind of instability develops. It has been shown in 
|K| that when the trap is not fully symmetric the first 
kind of instability exist for every direction of the angular 
velocity. 

We will show now, that the second kind of instability 
does not occur when the rotation vector is along one of 
the axes of the trap. Therefore, it can never happen in a 
two dimensional trap. Without any loss of generality, one 
may assume that the direction of rotation is along the z 
axis. In this degenerate case the characteristic polyno- 
mial Q factorizes because the motion in the z direction 
is not influenced by rotation, 

x(x-K) = 0. (8) 

The biquadratic polynomial in JHJ has real zeros if its 
discriminant A is positive and this is, indeed, the case 
since 

A = {2n^ + v^ + Vyf ~4{fi''~n^iv^ + Vy) + v^Vy) 
= sn^iv^ + Vy) -f (14 - Vyf > 0. (9) 



The second kind of instability always exists if the di- 
rection of rotation does not coincide with one of the axes 
of the harmonic trap. The tilting the axis of rotation 
away from the trap axis not only leads to the appearance 
of an additional region of instability but also, as shown 
in the next Section, leads to new effects when gravity is 
taken into account. 



IV. ROLE OF GRAVITY 

So far, we have neglected the presence of the gravita- 
tional field. It turns out that, owing to a possible reso- 
nance gravity may strongly influence the center of 
mass motion and lead to an instability. In the frame co- 
rotating with the trap the gravitational field is seen as a 
time dependent periodic force. This force, under appro- 
priate conditions, may fall into resonance with one of the 
characteristic frequencies of the trap and cause the lin- 
ear growth of the amplitude and the particles will escape 
from the trap. Of course, such a homogeneous external 
field will only infiuence the center of mass motion; it will 
not change the dynamics of the internal, relative motion. 
This is all true in the classical and in the quantum theory. 

The equations of motion for the center of mass in the 
presence of gravity acquire an additional term as com- 
pared to (0) 



dX(t) ^ , , / 
— ^ = MX(t) + 
dt ^ ' \mg{t) 



(10) 



In the rotating frame the gravitational acceleration can 
be written in the form |lll | 

g{t) = (gii + {gi_ + i{n x gi))e^"*) , (11) 

where g\\ and g±^ denote the parts of the gravitational 
acceleration parallel and orthogonal at t = to the di- 
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FIG. 4: K = 1, ^K, = 2, 14 = 3 = sm(^), Uy = 0, 
rtz — cos(^). Two horizontal lines enclose the lower region 
of instability. Both resonant frequencies lie in the lower region 
of stability. 

rection of rotation n. Thus, we have arrived at the clas- 
sic problem of a periodically forced harmonic oscillator 
in three dimensions. Resonances will occur when the 
frequency of the external force Vl coincides with one of 
the three characteristic frequencies uji of the oscillator. 
An interesting novelty here, as compared to the stan- 
dard forced oscillator, is that the characteristic frequen- 
cies themselves depend on $7. In order to find the con- 
dition for a resonance, one has to put w = in Q and 
solve for the resulting biquadratic equation 

DVL* + EVl'^ + F = Q, (12) 

where 

D = -2{Tt{V} -n-V-n), 

E = ^ Tr{V}^^Tr{y2} ^ Tr{T>}n-f-n- n- V^.^^ ^ 

F = -Det{y}. (13) 

As was shown in jllj | this equation always has two real so- 
lutions for ri^, the smaller root lies in the lower stability 
region. It turns out that we can choose the parameters 
of the trap in such a way that the larger root falls into 
the lower region of stability (Fig. , into the lower re- 
gion of instability (Fig. O , or into the higher region of 
stability (Fig. . These figures also illustrate the sim- 
plest graphical method to find the resonant frequencies. 
Since x = w^, in our figures [TO that depict the locations 
of the stability regions, the positions of resonances are 
found at the intersections of the parabola x — with 
the previously drawn curves. 

Full analytical description is very simple in the two- 
dimensional case, when the trap rotates around one of 
its principal axes. For definitness, let us assume that 
the is along the z axis whose orientation in space is 
still arbitrary. Without any loss of generality, one may 



FIG. 5: K = 1, 14 = 2, K = 3 = sin(f ), Uy = 0, 
Tiz — cos(-|). Higher resonant frequency lies in the lower 
region of instability. 
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FIG. 6: K = 1, = 2, V; = 3 = sin(^), Uy = 0, 
riz — cos(^). Resonant frequencies lie in two different regions 
of stability. 

assume that < Vy. In this case, Eq. (|12|l has two 
simple solutions. 

From Eq. iQ one can deduce that the lower region of 
stability lies below fi^ = 14, the region of instability is 
in the range Vx < < Vy, and the higher region of 
stability extends above Sl^ = Vy. One can also see that 
< Vx- Thus, the lower resonant frequency always falls 
into the first region of stability. The location of the higher 
resonant frequency varies. When Vz <Vx, then both fre- 
quencies lie in the lower region of stability. Vz > Vy then 
the higher resonant frequency falls into the higher region 
of stability. In the last case, when Vx < Vz < Vy the 
higher resonant frequency lies in the region of instabil- 
ity. This case in not interesting because in this region 
the amplitude grows exponentially anyway and the res- 
onance will not be noticeable. In the special case of the 
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rotation around one of the trap axis, the position of res- 
onant frequencies does not depend on the orientation of 
the rotation axis with respect to the vertical direction. 
In the general case, the calculations are much more com- 
plicated. 

V. CLASSICAL-QUANTUM 
CORRESPONDENCE 

A. Matrix Riccati equation 

Let us consider a state of a quantum-mechanical sys- 
tem with the Hamiltonian described by a Gaussian 
wave function of the form: 

* = C(t)exp(-ir.X(i).r), (15) 

where the matrix K is of course symmetric (its real part 
is positive definite) and C is a normalization constant, 

C{t) = \jBet{^{k{t))/V^}. (16) 

From the Schrodinger equation ihdt'i' — Ti.'^ with the 
Hamiltonian one obtains the evolution equation for 
the matrix K in the form of a matrix Riccati equation 
{h= l,m= 1) 

^ = -iK^ + ~ p, K]. (17) 
dt 

Mathematical theory of these equations is well devel- 
oped |JL2, 13J. In particular, the matrix Riccati equations 
find numerous applications in control theory [lj|- It is 
known that these nonlinear equations can be replaced 
by the linear ones. To this end we shall follow the proce- 
dure outlined by us previously [T^ in the same quantum- 
mechanical context and we shall search for solutions of 
Eq. lfT7|) in the form 

k{t) ^ -iN{t)D-\t), (18) 

which leads to 

- i — D'^ + iND-^—D-^ 
dt dt 

= iND-^NL)-^ +iV + iClNL)-^- mb-^n. (19) 

This matrix equation is satisfied when the matrices N 
and D obey the following linear equations 

^ = N~nb, (20a) 
dt 

— = -VD-QN. (20b) 
dt 

The linearization of the Riccati equation, in addition to 
being an effective mathematical tool, has also concep- 
tual advantages. Namely, it leads to a direct relationship 



between classical and quantum theory. Comparing H20|l 
with (0) , one can see that the columns of the matrices D 
and N satisfy the same equations as the classical position 
and momentum vectors, respectively. Therefore, from 
the knowledge of the classical motion, one may deter- 
mine the evolution of a Gaussian wave function. This re- 
lationship between classical and quantum mechanics has 
been noticed before by Arnaud [ig but only in the spe- 
cial case of the one-dimensional oscillator, when there is 
no nontrivial mode structure. Here, we extend it to three 
dimensions and we also give a prescription how to con- 
struct stationary quantum states from the classical mode 
vectors. 



B. Stationary quantum states 

In a stationary quantum state (|15|l the matrix N(t) 
must have the following form: 

N{t) = ikob{t), (21) 

where Kq does not depend on time. Inserting (|21() into 
(|20|l gives an obvious relation (cf. Eq. I17|l 

^ -ik^ + iV - [Cl, ko]. (22) 

However, we gained a new insight: stationary solutions 
of the Riccati equation H17|l can be constructed from the 
solutions of the set of linear equations (|2Jl of classical 
mechanics. 

In order to construct a stationary solution of H17I) for 
a harmonic oscillator in d dimensions we need d linearly 
independent eigenmodes of the classical equations of mo- 
tion. Such eigenmodes evolve in time according to the 
formulas Xi(t) = Xi exp(iWii). Let us construct now two 
dx d matrices x{t) and p{t) whose columns represent the 
positions and momenta, respectively. These matrices can 
be written in the form 

x{t) — xoe'Kp{iLut), p{t) — po exp{iLut) , (23) 

where w is a diagonal matrix whose elements are equal 
to the frequencies Wj of the chosen modes. Since the 
matrices x{t) andp(i) are solutions of Eqs. H20|) . one may 
identify them with the matrices N{t) and D{t). The 
matrix k — x{t)p{t)^^ will satisfy the matrix Riccati 
equation. Moreover, owing to the cancellation of the time 
dependent terms exp(zt2't), it will be a time- independent 
solution. 

In addition to the special stationary solution, there 
are also many nonstationary solutions of Eq. (|17|l . They 
describe pulsating Gaussians and they are the analogs of 
squeezed states, known from quantum optics. 

C. Two-dimensional harmonic oscillator 

In general, Eq. (|22|l for kg leads to a system of three 
quadratic equations that are hard to solve but one can 
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find the solutions quite easily in the simple case — for a 
the two-dimensional rotating harmonic trap (the Whit- 
taker problem). This explicitly soluble example will be 
used to explain the connection between the classical and 
quantum cases. 

We shall choose the coordinate axes along the direc- 
tions of the trap to make the potential matrix diago- 
nal V = diag{Vx, Vy,Vz)- The angular rotation vector 
is parallel to z axis of the potential ellipsoid. We shall 
parametrize the solutions of Eq. H22(l as follows 



Ko 



a 
ij P 
^/K 



(24) 



and obtain the following set of algebraic equations for 
a, l3, and 7 



2n-f + V, = Q, 



(25a) 
(25b) 
(25c) 



The solution of these equations proceeds in three simple 
steps. We first solve Eqs. (|25af) and (|25b|l for a and /3 
choosing the positive square roots to secure positivity of 
the real part of K 



(26a) 
(26b) 



Next we substitute these expressions into Eq. H25c|l to 
find 7 from the quadratic equation that results from 
squaring both sides of H25cl) . Only one solution of this 
quadratic equation satisfies (|25c|) and it can be written 
in the form 
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1 



(27) 



where k = zL^J {Vx — i^^)/ {Vy — il^) is a real parameter 
in both stability regions. In the first region the sign of 
K is positive, while in the second region it is negative. 
In the last step, we substitute the value (|27|l of 7 into 
Eqs. (|26() to obtain the final formulas for the parameters 
a and /3 

a = ^/Vx-n^ +4n'^/{l-l/K)^, (28a) 

p = ^Vy-n^ + Any{i-Ky, (28b) 

The matrix H24|l is symmetric, positive definite in all sta- 
bility regions and it satisfies the equation H22|) for sta- 
tionary states. 

In the three dimensional case a direct solution of the 
nonlinear equations for Kq is impossible since it requires 
solving high order polynomial equations. However, the 
method that employs a representation of K in the form 
(|18|l requires only the solution of linear equations. 



D. Construction in terms of the classical modes 

We shall show now that the general procedure de- 
scribed in IV Bl reproduces in two dimensions the solu- 
tion (|24|l . In the two-dimensional case there are four 
solutions corresponding to four characteristic frequencies 

2 Vx+Vy+ 2^(2 + ^{Vx-Vyy+Sn^iVx+Vy) 



Vx + Vy + 2n^ - ^{Vx - VyY + %^'^{Vx + Vy) 



/x{uj) 



N 



/ 2iLun 

Vx-Lu^- 



where iV is a normalization constant and is one of the 
four characteristic frequencies. 

From any two modes one may construct the following 
matrices: 



\y[uji) y{uj2)J 

jQ. ^ fpx{i^l) Px{uJ2)\ 

\Py{uJl) Py{uJ2)) ' 



(29a) 



(29b) 



and then construct = —iJ\f'D~^ . This matrix will 
satisfy the relation (|22|l for a stationary solution of the 
Schrodinger equation. There is still the condition that 
the real part of Kq is positive definite to guarantee that 
the solution is square integrable. One may check that the 
choice of two modes that differ only by the sign of the 
characteristic frequency does not give the right solution 
since in this case the real part of K has one eigenvalue 
zero. This leaves us with four possible choices of the 
signs of two different frequencies. When Vl lies in the 
first stability region {Q, < mm{Vx,Vy)) then both signs 
must be chosen positive. When lies in the second sta- 
bility region (f2 > max(T4;j VJ^)) then one must choose 
the opposite signs with the minus signs accompanying 
the lower frequency to. In the region of instability it is 
not possible to choose the signs in such a way that the 
real part of the matrix K is positive definite. The matrix 
K constructed in this way is unique (the normalization 
constants N cancel out) and coincides with the previously 
calculated matrix H24I) . The choice of signs is determined 
by the diagonal form of the classical Hamiltonian: 



Ti. — uJialai + ^20^0.21 



(30) 



where the complex amplitudes ai and a* satisfy the fol- 
lowing Poisson bracket relations: 



{a^, a*) = -i5ij {a,j, aj) = {a*, a*} = 0. 



(31) 



In the first stability region both terms in H3U|I are pos- 
itive while in the second region the term corresponding 
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to smaller value of w is negative. This is exactly the 
prescription used in constructing the solution 

Having found the wave functions, one may proceed to 
calculate the Wigner functions. They are especially use- 
ful for a comparison with classical physics. For Gaussian 
wave functions, the Wigner function is an exponential of 
a quadratic form of positions and momenta 



Wir,p) = Mexp{~^X-W-X) 



(32) 



For stationary states this quadratic form is a constant of 
motion and as such must be a combination of the expres- 
sions (IA2p . In the two-dimensional case it simply reduces 
to: 



W{r,p) = Mexp(fciCi + fcaCz), 



(33) 



where Ci and C2 are the constants of motion defined by 
[4|) and 



fci 

fc2 



2{(3Vy ~ aV^) 
al3{Vy 
2{a~ 



aP{Vy~V,)- 



(34) 
(35) 



The construction of the matrix K from the solutions 
of the classical equations of motions in three dimensions 
proceeds along similar lines. In this case the solution 
of the nonlinear equations H22I) by a direct method is 
more complicated, if possible at all. Our method of con- 
structing K from the eigenmodes is based on the fol- 
lowing unambiguous prescription, how to choose from 
the six modes the appropriate three modes. First, the 
three modes must have different absolute values of the 
frequency. The choice of the signs depends on the region 
of stability. In the first region, all signs are positive, in 
the second the smallest frequency must be taken with the 
negative sign, while in the third region it is the middle 
frequency. 

It is clear from our construction of the stationary Gaus- 
sian solution of the Schrodinger equation from classical 
modes shows that the classical and the quantum prob- 
lems always share the same regions of stability and in- 
stability. 



the Newton equations and the Schrodinger equation, the 
same regions of stability/instability are characteristic of 
both the classical and the quantum case. The dynamical 
properties of rotating traps are further enriched by the 
presence of gravity if the axis of rotation is not vertical. 
Gravity-induced resonances cause, at sharply defined ro- 
tation rates, an escape of particles from the trap. 

This research was supported by a grant from the Polish 
Committee for Scientific Research in the years 2004-06. 



APPENDIX A: CONSTANTS OF MOTION 

The Hamiltonian of every harmonic oscillator in three 
dimensions can be written as a sum of three independent 
terms 



H = ujia\ai + uj2a\a2 + 0^30303, 



(Al) 



where the coefficients are not necessarily positive. In 
the classical theory a's and a^'s are complex amplitudes 
while in the quantum case they become annihilation and 
creation operators. It is clear that each part of I|A1|1 is a 
constant of motion. These constants expressed in terms 
of p and r have the general form 



1 



1 



C = -prp + rVKp+ -rUr 



(A2) 



where the matrices T, W, U must have a rotationally in- 
variant representation in terms of ft and V. The time 
evolution of C, as determined by the equations of motion 
leads to the following equations for T, U, W 

= f ^ [n,f]+w + 
O^U ^ [n,Lf]-WV -VW^ (A3) 
o^w = [Cl,W] + U -vf 

The number of independent quadratic constants of mo- 
tion is the same as the number of dimensions of the har- 
monic oscillator. They can be chosen as the bilinear com- 
binations of the mode amplitudes a^Ofc. 
In two dimensions we have 



VI. CONCLUSIONS 

We have presented a full analysis of the dynamical 
properties of an anisotropic, rotating harmonic trap. The 
main result of this analysis is the discovery of three dis- 
tinct regions of stability and two regions of instability. 
In the generic case, by increasing the rate of rotation we 
visit successively all these regions. We have also shown 
that, owing to a direct link between the solutions of 



Ci 



Co = 



f = 


t 






w = 


n 






ii = 


V 






f ^ 


V 






w = 


nv- 


f 2T>f2- 




u = 




- 





(A4a) 
(A4b) 



In three dimensions there are three independent solu- 
tions of Eqs. 



J 
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T = t 

Ci = { U = V (A5a) 

V ^ h 

f = V-iVL^ 

u ^ - vn"^ -Ci'^v ~ nvn (A5b) 

^ " -wn^v - evCivn - eCivtiv ^ bn^v^ + i3Tt{vCi^}v ^ 



The first constant l|A5a|l gives the total energy (Hamil- 
tonian) of the system . The remaining constants are 
higher order expressions in V and Q and they have no 



direct interpretation. The energy of each mode of the 
oscillator is a linear combination of these constants. 
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